Superscars in billiards — A model for doorway states in quantum spectra 
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In a unifying way, the doorway mechanism explains spectral properties in a rich variety of open 
mesoscopic quantum systems, ranging from atoms to nuclei. A distinct state and a background 
of other states couple to each other which sensitively affects the strength function. The recently 
measured superscars in the barrier billiard provide an ideal model for an in-depth investigation of 
this mechanism. We introduce two new statistical observables, the full distribution of the maximum 
coupling coefficient to the doorway and directed spatial correlators. Using Random Matrix Theory 
and random plane waves, we obtain a consistent understanding of the experimental data. 
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Strength function phenomena [l[ in open mesoscopic 
quantum systems are a central object of study in atomic 
and molecular physics as well as in atomic clusters, quan- 
tum dots, and in nuclear physics . Often there is a some- 
how "distinct" and "simple" excitation whose amplitude 
is spread over many "complicated" states. The distinct 
state thus acts as a "doorway" to the (usually chaotic) 
background of the complicated states [l|, 0|. Prime ex- 
amples are Isobaric Analog States and multipole Giant 
Resonances (GR) in nuclear physics. The strength func- 
tion is typically of Breit-Wigner (B W) shape with a char- 
acteristic spreading width 0, 13, S IB] • For examples from 
molecules and metal clusters, see 

Quantum billiards can be realized experimentally by 
flat microwave resonators [T3|. To study the doorway 
mechanism in detail, we use a microwave billiard of rect- 
angular shape with a thin barrier inside, see Fig. [TJ The 
electric field strength distribution corresponding to the 
quantum wave function is reconstructed from the mea- 
sured intensities. Certain wave functions of this pseu- 
dointegrable billiard possess unique structures called "su- 
perscars" [ill . [l2| . These are scarring wave functions re- 
lated to families of neutrally stable classical periodic or- 
bits. Four examples of measured superscars are shown 
in Fig. [TJ Unlike ordinary scars [13|, which are local- 
ized around a single unstable periodic orbit, they do not 
disappear at large quantum numbers. They are embed- 
ded into, but clearly distinct from, a large number of 
nonscarred wave functions. We will demonstrate that 
the superscars act as doorways to the background of the 
nonscarred wave functions. Our perfect control over the 
experimental observables allows us an in-depth study of 
the doorway mechanism which can presently not be ac- 
complished in traditional quantum systems. 

First, we briefly compile the necessary information on 
measured superscars in the barrier billiard. Second, we 
introduce the maximal coupling coefficient as a new ob- 
servable and use Random Matrix Theory (RMT) 
to model its distribution. Third, we introduce directed 
spatial correlators as another new observable and model 
them by extending Berry's random wave ansatz (l5| . 




FIG. 1: Examples for measured superscars induced by the 
barrier and concentrated along the indicated classical periodic 
orbits (dashed lines). They are members of four different 
families. Top row: horizontal Bouncing Ball BB and Inverted 
V superscars; bottom row: Diamond D and W superscars. 
The gray level indicates the value of the wave function (black: 
highest positive and white: most negative value; see [12l|). 

As Fig. Q] shows, the superscars with a clear wave func- 
tion structure relate to particular classical periodic or- 
bit. The superscars form a family which lives within 
an infinitely long Periodic Orbit Channel (POC). Due 
to diffraction on the tip of the barrier, the amplitude of 
the scarred wave function tends to zero along the POC 
boundary. Thus, the superscarred wave function can 
be approximated by a constructed superscar state, de- 
fined as an eigenfunction ^m]n(r) in the infinitely long 
POC [H[I3. Here Fe{BB,V, D,W} stands for the su- 
perscar families, as defined in Fig. [TJ and (m, n) are the 
numbers of wave maxima along and perpendicular to the 
POC. A measured state ^ f(r) at (rescaled) frequency / 
in the barrier billiard has an overlap, 

Cm,n = (*£?„|*/) , (1) 

with the constructed superscars [lj], E2]- As an exam- 
ple, the distribution of the overlaps with a constructed V 
superscar state with quantum numbers (m, n) = (45, 1) 
is depicted in Fig. [2j The superscar strength spreads 
into a neighboring nonscarring background states follow- 
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ing a BW shape with the main strength concentrated in 
a few states. This nicely confirms our doorway interpre- 
tation. For comparison, a nuclear GR doorway is also 
plotted in Fig. [2 Here, the number of background states 
— reflected in the fluctuations around its BW shape — 
is much larger than in the barrier billiard. According to 
the Brink- Axel hypothesis [lB,[l3| a GR excitation builds 
upon every nuclear state. Similarly, a superscar doorway 
state exists for each value of F, m and n. 

We now set up a random matrix model in the spirit 
of models in nuclear physics [3, 0] . The total Hamilto- 
nian reads H = H s + Hb + V . Here, H s and Hb describe 
doorway states and background states, respectively, and 
V couples the two classes of states. The eigenequations 
for the uncoupled Hamiltonians are H s \s) = e s \s) and 
Hb\b) = eb\b). For the matrix elements of the interac- 
tion, we make the assumptions (s|V|s') = (&|V'|6 / ) = 
and (6|y|s) = Vb s for any s, s', b, b' . We interpret the 
constructed superscars, 

*$,n(f) for a given family F but 
with different (m, n) as the doorway states s. Due to the 
interaction V the doorway state is not an eigenstate of 
the Hamiltonian H. We assume that the interaction ma- 
trix elements, Vb s — v s b, are Gaussian distributed random 
variables with zero mean and variance v 2 . Importantly, 
the parameter governing the physics isv/d, where d is the 
mean level spacing of the background states [l|, ll4( . Since 
only a few states carry superscar strength with given val- 
ues of (m, n), and superscar states with different (m, n) 
are assumed to not mix, it is sufficient to consider only 
one superscar state, s, coupled to N background states, 
b, where N is large. To resemble the experiment, we in- 
clude iV=294 background states. As the barrier billiard 
is pseudointegrable, the spacings between the eigenstates 
are semi-Poisson distributed We thus generate such 
an ensemble of N + l states. The doorway state is chosen 
as the middle state and interacts with the surrounding TV 
states. For each realization, energies and wave functions 
are numerically obtained, and the mixture of the super- 
scar with the surrounding states is calculated. We then 
extract v/d for each superscar family. The full problem 
H\n) — E n \n) is solved by the exact implicit equation 

N 2 

£» = es-£-^, (2) 
b=i eb ^ n 

and the wave functions are given by 

I»>=*(|*>- ■ ( 3 ) 
The superscar coupling of each eigenstate is therefore 

/ * 2 \- 1/2 

C * (n) =l 1+ £(^M ' (4) 

The superscar strength c? over the different eigenstates 
\n) is BW distributed [l| with spreading width T^- = 
2Trv 2 /d,i.e.r l /d=2Tr(v/d) 2 . 
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FIG. 2: Doorway strength functions. Left: overlap between 
constructed superscar state of the V family with m = 45, n = 
1 and the measured wave functions versus rescaled frequency. 
The solid curve is a BW function. The inset shows the overlap 
on a logarithmic scale over a large frequency interval. Right: 
spectrum of the (p,p) reaction at 200 MeV on 120 Sn in the 
region of the Isoscalar Giant Quadrupole Resonance (ISGQR) 
as an example for a nuclear doorway state The solid curve 
is a BW function fitted to the data. 

As the superscar strength is distributed over a small 
number of states only (see Fig. [2| , v/d is smaller than 
or of the order of unity. The fit with a BW distribu- 
tion shows a rather large variation of the fitted shape as 
well as the width of the distribution over the ensemble of 
observed superscars in a superscar family F. Hence, the 
width is not a well-suited measure to determine v/d. 
We thus consider the state with the largest coupling 

c ma X = max (<4, n ) (5) 

for a given superscar which is directly obtained from ex- 
periment. Since a rather small number of states carry 
strength from the doorway state (i.e. the constructed 
superscar), the peak of the fitted BW shaped strength 
function usually deviates from the measured largest su- 
perscar strength: The discretely measured state does not 
appear exactly at the peak. We may, however, directly 
compare the maximal measured value to the correspond- 
ing calculated value, ma,x(c 2 (n)), where c s is obtained 
from Eq. (@]) . Not only the average value of c^ ax can be 
studied but also its higher moments. We study the full 
distribution of these maximal couplings for a superscar 
family F, which, as far as we know, has never been con- 
sidered before. In Ref. the first two moments of the 
c max distribution were studied, but with assumptions not 
valid in our context. The shape of the c^ lax distribution 
strongly depends on the interaction strength, v/d, and 
it is a particularly sensitive measure for small values of 
v/d, i.e. of the order one or smaller. 

In Fig. [3] we show measured distributions of c^ ax with 
the best fit curves of the RMT model for each superscar 
family F. The fit gives the following values for the interac- 
tion strength: for the BB superscar v/d — 0.45, for the V 
superscar v/d — 0.35, for the D superscar v/d — 0.3, and 
for the W superscar v/d = 0.55. The coupling strengths 
are small and thus our ansatz for a BW shape for the 
doorway strength function (Fig. [2]) is in accordance with 
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FIG. 3: Left: c^ax distributions of measured superscars (his- 
togram), and the fit of the RMT model predictions (solid line). 
Right: normalized distributions of superscar strength spread 
over all states on a logarithmic scale. Experimental distribu- 
tions (dots) are compared with the RMT model predictions. 



earlier findings |20|. The V, D, W superscar families con- 
tain 16, 25 and 22 measured members, respectively, while 
the BB superscar family contains only 9. The fit in this 
latter case has thus higher uncertainty. The averaged 
measured and calculated Cj 2 nax values are listed in Tab.U 

Another observable is the distribution of the super- 
scar couplings over all eigenstates. The strength of each 
constructed superscar is measured (and calculated) over 
all 294 states, where the major part of the strength is 
concentrated in a few states only. Figure [3] shows mea- 
sured distributions compared to calculations for different 
interaction strengths obtained from the fit to the c^ ax 
distributions. Once more, we clearly see that the model 
reproduces the experimental distributions for all super- 
scar families well except in the case of the BB superscar 
family because of the small number of superscars. 

We now turn to the spatial correlations of the wave 
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0.9 ±0.1 


1.3 


V 
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0.74 ± 0.03 


0.72 0.69 


0.9 ±0.1 
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0.54 ± 0.03 


0.51 0.49 


1.0 ±0.1 


1.9 



TABLE I: Experimental results with standard errors of the 
mean versus results from RMT model and directed correlators 
(Corr) for averaged c^ax values and spreading width . 



functions. Berry [15|| introduced the correlator 

(MR-mr k (R+r/2)) 



C{kr) 



(\MR)\ 2 ) 



(6) 



of the wave functions ipk(r) where the average is per- 
formed isotropically over all vectors R and, for fixed mod- 
uli of wave vector k and r, over all directions of the vector 
r. In our context, all wave functions are real and no com- 
plex conjugation is needed in the definition ((6]). Berry 
argued that the spatial correlations of a wave function 
in an ergodic system should be indistinguishable from 
those of superimposed plane waves. In two dimensions 
this yields the universal prediction C(kr) = Jo (At), if 
possible boundary effects are ignored. Here Jo is the 
Bessel function of order zero. Indeed, this behavior was 
confirmed in numerous systems [l(J E3, E3, H2] . 




FIG. 4: The wave function correlators. The Jo(kr) predic- 
tion is always given as dashed line. The top row shows the 
correlation function of the constructed V superscar state as 
solid lines: the isotropic C(kr) as well as the directed C ± (kr) 
and C"(fcr). In the middle row, the same observables are de- 
picted as solid lines for the averages over all experimental wave 
functions in the barrier billiard. In the bottom row, the cor- 
relators averaged over all observed V superscars are displayed 
as solid lines and the correlators resulting from Eq. J?} with 
Cmax =0.69 are shown as filled circles. 

The superscars, however, clearly have non-ergodic fea- 
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tures. To analyze their correlations we define new, es- 
pecially tailored observables to which we refer as di- 
rected correlators. Instead of averaging isotropically as 
for C(fcr), we now carry out the averages either only 
across or only along the channel in which the superscar 
exists similar to [23|. We thereby obtain the correla- 
tors C ± (kr) and C"(kr), respectively. In the top row 
of Fig. H] the three correlators of a constructed V super- 
scar ^m^n(r) are depicted. While the isotropic correlator 
C{kr) follows the Jo(fcr) prediction up to a certain scale, 
the directed correlators strongly deviate from it. The re- 
sults for C ± (kr) and C"(fcr) show that the constructed 
superscar fills the channel and moves through it as a sine 
wave, see also Fig. [H This information about the form of 
the waves, however, is washed out when averaging over all 
wave functions in the billiard. As the middle row of Fig. [4] 
shows, each of the three correlators worked out for all 
measured wave functions coincides with the Jo(fcr) pre- 
diction for chaotic systems. Hence, we may use Berry's 
random wave approach even though our billiard system is 
pseudointegrable. Importantly we only use the two-point 
correlations and only go up to kr — 8. 

We now use these observations to extract information 
about the superscar couplings from the measured corre- 
lators. Correlators averaged over all experimentally ob- 
served V superscars are displayed in the bottom row of 
Fig. [H They are similar to, but slightly different from 
those for the constructed V superscars in the top row. 
The difference is due to the leaking of the superscar out 
of the channel or, in the language of the doorway de- 
scription, due to the coupling of the background states 
to the superscar. We thus model the measured super- 
scars ^p(r) for family F as a linear combination of a 

constructed superscar ^m] n {r), which only contributes 
in the channel, and a state Xfc(r) which is ergodically 
distributed everywhere in the billiard, 

* ( jP(0 = <W*£?n(0 + V^&Zxktf ■ (7) 



This ansatz is fully consistent with the RMT model set 
up above and extends it by also modeling the spatial de- 
pendence. The states describing the background should, 
first, have Jo(kr) correlations and, second, be orthogonal 

to ^m] n (r). Thus, we choose the "scarless" plane waves 



Xk(r) = , (8) 

with standard plane waves Xfc(j)- The superscar con- 
tribution in the plane waves is small (but not negligi- 
ble); the distribution of the overlaps (^m] n \xk) has a 
standard deviation of 0.13. We convinced ourselves that 
the correlator of the Xfc(r) follows the Jo(kr) prediction 
very closely. We work out the three correlators for the 
model (0. They depend on c max which is, just as in the 
RMT model above, the coupling to the superscar door- 
way. By fitting to the measured superscar families we 
determine the couplings c max . The fits for the V super- 
scar are shown in Fig. |H The resulting (c^ ax ) values in 
Tab. U are close to those obtained from the RMT model. 
This is a nice mutual confirmation. For comparison, we 
also give the resulting V s - values in Tab. [H Obviously, 
our new observables are more appropriate. This is born 
out in the large standard deviation of the 1^ distribution 
which is, e.g. for the W superscar family, 0.8. 

We conclude that our doorway interpretation yields 
a thorough understanding of the experimental findings. 
Our two new observables give deeper insight into the sta- 
tistical features of the doorway mechanism as such, and 
it is encouraging to see how well the two analyses agree. 
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